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Abstract 

We examine the effect of adding pion exchange between quarks and pion self-energy corrections 
to the Los Alamos Relativistic Quark model with a short-distance cutoff of the Bethe form. The 
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under this change in the sense that significant changes to the model parameters are not required. 
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I. INTRODUCTION 



The Los Alamos Relativistic Quark model (LARQ) has been applied to the ground states 
of octet and decuplet baryons and to small nuclei (^He and ^He jsl) with good 
numerical results. The motivation is the conjecture that, just as distorted atomic electron 
wave functions can replace bond energies, electron affinities, and electronegativity, etc., and 
so provide chemistry with a quantum mechanical description of the binding of atoms into 
molecules, so too, a corresponding distortion of quark wave functions can provide an accurate 
description of the binding of nucleons into nuclei, and replace previous descriptions in terms 
of internucleon potentials or exchanges of mesons. The approach also implicitly provides a 
description of "off-shell" nucleons, or equivalently, of the structure of a nucleon in a nuclear 
medium. 

For isolated baryons, the model has determined current quark wave functions by solving 
the Dirac equation with a linear scalar confining potential, although in principle, any other 
(relativistic) method for the calculation could be employed. In addition, the spatial currents 
of single gluon exchange are evaluated with these wave functions to calculate the color- 
magnetic spin-spin (CMSS) interaction to the total state energy in its rest frame, which are 
essential to distinguishing the nucleon from the A-baryon. A Gaussian spatial propagator 
with a range of order 1 fm is used for these gluon exchanges. 

For nuclei, a geometrically complex scalar mean field potential is introduced, consisting 
of an appropriate array of the scalar potentials truncated at the mid-planes between the 
array points. Variational wave functions are then employed, which consist of amplitude- 
weighted sums of the quark wave function solutions for the isolated potential, centered on 
each of the different array locations. The quark amplitude is assumed to be dominated by 
the contribution from a single array location (potential well) with the relative amplitude 
of the quark wave functions from the other wells taken as a variational parameter. Three 
quark wave functions are color-spin-isopin correlated to nucleon quantum numbers for each 
(dominant) well. The distributed coherent quark amplitudes are interpreted as a delocal- 
ization of the quarks from each nucleon. A final overall antisymmetrization is implemented 
to satisfy the Pauli principle. As well as the relative (delocalization) amplitude, the scale 
for the separation of the array points is used as an additional variational parameter to find 
the minimum energy configuration (in a body-fixed frame). A variational estimate of the 



single-body energy is obtained as the square root of the matrix element of the squared Dirac 
Hamiltonian for the multi-well potential with these wave functions and the contribution of 
single gluon exchange to the total energy of the state is evaluated as discussed above. 

The effect of delocalization is as one would expect on general quantum mechanical 
grounds: the kinetic energy of the quark is reduced and so likewise the overall energy of the 
state. It is of little consequence both when the potential wells are very close together, and 
also when they are far apart. In the former case as there is little additional volume for the 
delocalization to occur, and in the latter the amplitude suppression from tunneling 

through the confining barrier formed by the truncated potentials becomes severe. Even at 
an intermediate well separation, however, delocalization does not proceed unimpeded. This 
is because the mass of a nucleon is reduced by the CMSS interaction (while it is increased for 
the A-baryon). As this interaction decreases with decreasing quark density, the net effect 
of delocalization here is an increase in the energy of the bound nucleon state corresponding 
to a repulsive interaction. These two opposing trends reach an equilibrium in the nuclear 
calculation and produce a modest distortion of the quark wave function for a bound nucleon 
relative to that for a free nucleon. 

The binding energy found for ^He is 4 ± 2MeV and for ^He is 20 ± 6MeV. The (uniform) 
well separations at these energy variational minima are approximately 1.8 and 1.75 fm 
respectively. The relative amplitude parameters are ~ 0.10 and ^ 0.14, corresponding to 
only a few percent difference between the probability distribution for a quark in the nuclear 
wave function compared with that for an isolated nucleon. 

One thus has a dynamical confirmation that, despite the huge strength of the strong 
interactions, something quantitatively identifiable as nucleons are the quasi-particles that 
build the nucleus, as has been assumed since the discovery of the neutron. This is analogous 
to the construction of molecules in terms of identifiable atoms despite the distortion of the 
atomic electron structure from that in isolated atoms. 

When the nuclear quark wave function for ^He is applied to calculation of the valence 
quark structure function for that nucleus, the resulting EMC effect \^ (deviation from one of 
the ratio of the deep inelastic lepton scattering response function of the nucleus to that for A 
nucleons) is found to be consistent [sl with the analysis by Kumano 6] of the A-dependence 
of the EMC effect. This occurs with no free parameters available for adjustment to match 
the data. 



A non-relativistic version of this model, the Quark Delocahzation, Color Screening Model 
(QDCSM), employs harmonic oscillator confining potentials that are screened [3] instead of 
truncated. It reproduces known baryon-baryon potentials [8] as far out as the intermediate 
range (~1.5 fm). This includes the well-known intermediate-range nucleon-nucleon attrac- 
tion, and reproduces the experimental ^Sq scattering phase shift j^. At larger distances, 
however, the interaction potential decreases at a Gaussian rate rather than exponentially, 
due to the Gaussian form of the quark wave functions. Since the scale for the Gaussian 
fall-off is ~ 1 fm, the model does not describe the deuteron at all well. An extension to this 
.node, was ,„,p.e,„e.ted Q by add,.g pion exchange between ,ua* w.th a *o.t-d,stance 
cutoff. The rationale for the cutoff is that all meson exchanges at short distances are already 
taken into account by delocahzation of the quark wave functions. With a cutoff of order 
0.6 — 1.0 fm, a deuteron of approximately normal size is obtained with a value for the ratio 
of D- to S-wave components of 4%, comparable to the normally reported 5%. Improvements 
to the ^Si and ^Pi phase shifts are also found. 

We therefore consider here the effect of adding pion exchange between quarks to the 
relativistic LARQ model, also with a short-distance cutoff, of the Bethe form This 
effectively gives quarks a form factor. For consistency with the (broken) chiral symmetry of 
Quantum ChromoDynamics (QCD), we must also include the change in the self-energy due 
to pion-loop corrections to the quark propagator. The form factor amounts to converting (for 
this purpose) the current quarks of the LARQ to constituent quarks with a finite structural 
size: the energy takes the place of the mass although the wave functions remain unaltered 
from those of the massless current quarks. 

The outline of the paper is as follows: In Sec. II we calculate the contribution of pion 
exchange between different quarks in the baryon wave function. In Sec. Ill, for consistency 
with chiral theory, we calculate the contribution of pion loops to the quark self-energy in 
three approximations: 1) no excitation of the bound-state quark to a higher mode (as in 
Ref. 121]), and 2) two variations with a free quark propagator in the intermediate state. The 
first allows comparison with previous work. As none of these results is large, we conclude that 
the size of this contribution is not significant. Finally, in Sec. IV, we present a discussion, 
including the relation to chiral perturbation theory, and our conclusions. 



II. PION-EXCHANGE BETWEEN DIFFERENT QUARKS 



The contribution to the self energy of the nucleon (or the A resonance) due to the 
exchange of a pion between any two of its three valence quarks, labelled "1" and "2" , is ^] 

i^S^'Hw) = -9',. I dh, I rfV2^^|^0(r - reut)(^i^75r,^i)(^2^75rfc^2) • (1) 

Note the lack of a "crossed" term corresponding to 1 — 2 and 2—^1. This is because 
the quarks in a baryon are distinguishable, necessarily having different colors. (We suppress 
here, and in the following, the antisymmetric color wave function.) We shall discuss the value 
of the coupling constant g^^^ in detail below. The pion propagator depends on r = |ri — r2|, 
the distance between the two quarks. The S-wave quark wave function ipi has the form 



Air I 



, (2) 



^ -icr ■ fi ipbin] 

with spin and isospin functions x ^"^^ ^- A similar expression holds for tfj2. 

As the LARQ models the contribution of all meson-exchange effects at short distances, 
we are only interested in finding the contribution from pion exchange at large distances. 
Equation ([T]) contains a G-function to assure that the six-dimensional integral there is to be 
evaluated for r > r^ut ~ 1 f^i- For example, the most recent Nijmegen jlj] fits to the A^A^ 
scattering data cut the unambiguous pion-exchange contribution to the analysis at Vcut = 
1.4 fm. 



A. Quark Wavefunctions 



The radial wave functions il!a,b in Eq- (E]) were obtainec. 
massless quarks in a linear scalar potential of the form 

V{r) = K^{r — tq) , 



w solving the Dirac equation for 

(3) 



where the negative offset, —k^vo, effectively gives these massless quarks a constituent mass. 
(It also affords inclusion of a rough representation of the effect of the short-distance color 
Coulomb interaction between quarks.) For the values used in Goldman, Maltman, Stephen- 

n 

son and Schmidt jSj] (referred to henceforth as GMSS), = 0.9 GeV/fm and tq = 0.57 



fm, the bound quark has a total energy of 361.6 MeV . The upper and lower radial wave 
functions, ipai'f') and ipbif), are chosen real in this description. 

The radial functions ipa,b{i^) were found in GMSS by numerically solving the Dirac equa- 
tion with the linear scalar potential of Eq. ([3]). In this form, they have dimensions fm~^/^ 
and are normalized so that 

POO 

Sr ^Ij\v)^{y) = Air / r'^dr [ipUr) + '^'^{r)] = 1 (4) 







Following Critchfield [15|], we can express the ipa,b{^) in terms of dimensionless, unnormalized 
functions 0a,fc(/tr), where k, = 2.1357 fm~^ (the GMSS value), as 

^,,,(r) = N «:3/20,,,(/€r) . (5) 

Explicitly separating out the natural scale of the problem, k, allows us to define the nor- 
malization factor, A^, to be dimensionless also. 

The numerical solutions for (f)a,b{x), x = K,r, have been well fitted (in GMSS's Appendix 
B) as a sum of Gaussians, 

n 

4>a{x) = ^ Ofc^ exp(-/ifc^a;V2) , (6) 
fci=i 

n 

(f)b{x) = xY^ 6fc2 exp(-/ife2a;V2) , (7) 

fc2 = l 

where we have defined fik = /ife + 1. The (dimensionless) coefficients a^, bk, and fik are 
tabulated in Appendix B of GMSS for n = 6, 12, and 18 Gaussian terms. For our purposes, 
the table for six terms provides sufficient accuracy. With Eqs. (jH]) and ([7j) for (f)a,b, we can 
calculate the normalization constant, finding 

roo 

Ar-2 = 4vr / x^dx [(f)l{x) + (pHx)] = 170.37 = (13.05)^ . (8) 
Jo 

B. The Interaction Energy 

Let us now calculate the interaction energy Bj^^'^'' (rent) between these two quarks. With 
Eq. ([2]), the two wave-function factors we need are 

^2h5rk^2 = '2^IJa{r2)Mr2) r2,/3 (xl^^ (T/S XS2J ^h,) ■ (10) 



Thus, 



B^^^ = -2 (^) ixl^a^ XsJ {xl,cr, X.J i^lr 6J ■ (4,t 6J 4, , (11) 



where 



lapircnt) = / d^ri / (i^r2 fi,„f2,/3 V^a(ri)V^6(ri) <S>{r - rent) ipa{r2)ipb{r2) 



-fix 



= N K J d xi J d X2 xi,aX2,/3 (l)a{xi)(f)bixi) ——Q{x - x^^t) 0a(a;2)06(a;2) ■ (12) 

Here fi = rriT,/ k, Xi = Kri, X2 = Kr2, a; = |xi — X2I, and Xcut = ^'^cut- 

Since all vectors in Eq. (fT2l) are integrated out, the tensor lap reduces to 

Iaf3{r cnt) = Sap /(rcut) , (13) 

where 

H^cut) = J d^Xi J d^X2 Xl ■ X2 4>a{Xi)(f)b{Xi) ——Q{x - X^nt) 0a (3^2)06 (2:2) • (14) 



Physically, this simplification happens because the pion-exchange interaction between two 
spin-i particles in a relative S'-wave can only involve the central scalar, but not the tensor, 
contribution. 



C. Evaluating the Integral /(rcut) 



To evaluate the integral of Eq. (IT4l) . we will do the X2 integration first, rewriting 

Hrcut) = / d^xi 0a(a;i)0fe(a;i) J(xi,Xcut) , (15) 



where 

-fix 



d^X2 0a (3^2) 06 (3^2) Xl ■ X2 Q{x - X^nt) ■ (16) 

X 

To proceed, let the integration variable for J(xi,a;cut) be x = xi — X2, with x running from 
Xcut to oo. As Eq. f|T6l) has no vectorial dependence remaining after the integrations. 



J(xi ) = J{xi 



= 27r/ x^dx dcose ^ x^ e 0^(|x, -x|)0b(|xi-x|) , (17) 

J Xcut J — 1 |Xx — X| X 

with cos ^ = Xl ■ X. Defining 

Cfcifc2 = (/^fci +/^fc2)/2 (18) 



and using Eqs. and ([7j), 



6„(|xi - x|)05(|xi - x|) = |xi - x| ^ ak,h,e-'"'i''^\^'-^^" . (19) 



The factor of |xi — x| here cancels the denominator in Eq. f|T71) . Expanding |xi — xp as 
x1 + x"^ — 2xiX cos 9, we can write 

J(Xi,Xcut) = X! ^k-fik2 Jkik2{^l7 ^cut) , (20) 
fcl,fc2 

where 

JkM^u^cnt) = r dx [' dcosO x{xi -xcos9) e-/^--^i'=2(-?+-'-2-i-™se) _ ^^l) 
The double integral in Eq. (12T|) can be done analytically, yielding 

<^fcife2(^l? ^cut) = ■^"f [^'^I'^a (-^l? -^cut) ~ -BA;ifc2(a;i, Xcut)] 5 (22) 



where 



'^felfc2 

and 

1/2 



5..fc2(^i,^cut) = (^] e+'^'/^^'=i'=^{(l + /iXi)e-'^-^erfc[4tl2] 

-(l-/.xi)e+'^^^erfc[tH]}. (24) 



16j.) In the equation for Bkj^k2 "we have 



(For definitions of the error functions, see Ref. 
defined 

4tfe2(^l' ^cut) = y/ci^{Xcnt T Xi) + fl/ {2y/c^) . (25) 

One might worry that the integration over xi that follows in Eq. (|T5l) would blow up because 
of the e"'"'^^^ in the second term of Bkik2- However, this exponential growth will be damped 
by the Gaussians in the 0(j(xi)0b(xi) factor. 

We have at this point reduced the six-fold integration in Eq. (ITT]) to 36 simple, compact 
integrations over xi, which we can easily (and quickly) evaluate numerically. 

I{r^nt) = — ^ — J2 (^k,bk2 4ifc2(^cut) , (26) 

fcl,fc2 = l 

where 

roo 

hikiircut) = / dxi (f)a{xi)(f)b{xi) [Akik2{xi,Xcut) - Bk^k2{xi,Xcut)] ■ (27) 
Jo 



We have done the numerical integrations using Mathematica [17|, which has a built-in 
erfc(x) function. For rcut = 1-0 fm, for example, / = —0.487 MeV. The dependence of 
-^('"cut) on Tcut is displayed in Fig. [TJ Note that this function is repulsive for rcut < 0.6 fm 
and attractive after that. This can be understood as due to the interplay between the cosine 
from fi ■ 02 and the cutoff function 0(r — rcut)- When rcut is small, the angle between fi 
and T2 can also be small enough so the cosine is positive. This combined with the resulting 
smaller r means that the positive pion propagator factor can give larger contributions to the 
integrand. However, when rcut gets sufficiently large, the two vectors will tend to contribute 
when they point in opposite directions. This makes the cosine negative, which makes the 
integrand negative for rcut > 0.6 fm, even though the pion propagator factor is diminished 
in size. 

As a check on our algebra, we have also calculated the triple integral in Eq. ( ITSi) for a 
few cases of rcut numerically. This takes, of course, considerably more computing time and 
the calculation is slow to converge. However, the numerical results are in good agreement 
with the results shown in Fig. [H 



D. Spin-Isospin Dependence of the Interaction Energy 

Altogether, contribution to the nucleon or A mass from the interaction energy of cut-off 
pion exchange between quarks 1 and 2 is 




with the spin-isospin factor 

x^'''^ = xsj ■ {xl,o- XsJ (el,,T 6tJ ■ (4,^ 6J • (29) 

Summing up over all pairs of the three quarks, the total vr-exchange interaction energy 
contribution, after taking into account the coupling of the quark spins and isospins to total 
S and T, is 

Bit = -4 ) Hrcut) XsT , (30) 

where 

XsT = 5] < S, Ms\(T^'^ ■ (T^^^\S, Ms >< T, Mt\t^'^ ■ t^^^\T, Mt > . (31) 

i<j 



The spin-isospin sum Xst can be evaluated in several ways. The most direct (most 
simplistic) is to write the operators Tj ■ Tj and cr^ ■ cTj in terms of their raising, lowering, and 
the 2;-components. Then apply these expanded operators, in turn, for each combination of 
i < j on the three-quark wave functions for the nucleon and A. This proceeds quickly for 
the A state of highest weight, | A"^"*", +3/2 >, and gives, after summing, the value 3. For the 
proton with spin-up, \p, +1/2 >, however, it is a rather tedious calculation which eventually 
yields the value 15. 

A more group-theoretic way of arriving at these numbers is given in an Appendix. How- 
ever, a simple, heuristic derivation of Xst goes as follows. First note that, for any pair of 
quarks, 

< (Ti ■ (Tj >= 4 < Si ■ Sj >= 2^(^ + 1) - 3 , (32) 

which is equal to +1 when the two quarks are in a triplet state and —3 when in a singlet. 
Likewise, Ti-tj is +1 and —3 in iso-triplet and iso-singlet states. As mentioned earlier, each 
S'-wave quark pair must be symmetric in the spin-isospin space. 

For the A, therefore, to get the total spin and isospin to both be 3/2, each quark pair 
must have S = T = 1. Thus 

< (Ti ■ O-j >< Ti ■ Tj >= ( + 1)( + 1) = 1 (33) 

and summing over the three pairs then yields 3 for the value of A's Xst- 

For the nucleon, however, each quark pair can be in either anS' = T = loranS' = T = 
state, and it is equally likely to be one or the other. Thus 

< CT, ■ (T, X T, ■ T, >= i(+l)(+l) + ^(-3)(-3) = 5 . (34) 

Now summing over the three pairs yields 15 for the value of nucleon's Xst- 

Note that, between two quarks, the spin-isospin factor is either +1 (S* = T = 1) or +9 
(S* = T = 0). Thus, for rcut > 0.7 fm (see Fig. [T]), Eq. (l30ll shows that the net interaction 
energy, Bj^l^'^^ , is positive. This result is to be contrasted with the case of two nucleons in an 
S-wave. In that case there is no color quantum number to provide for antisymmetrization. 
Thus antisymmetry of the two fermions involved requires that S = and T = 1 or = 1 
and T = (deuteron), and hence the interaction energy from the pion exchange has the 
opposite sign: (—3) x (+1) = (+1) x (—3) = —3. Thus, while this pion-exchange interaction 
is attractive in both the deuteron and threshold state, it is repulsive between quarks in 



both the nucleoli and the A-baryon, as one expects for a self-energy contribution to a stable 
particle. 



E. Evaluation of the Full Exchange Energy i?*°^(rcut) 

Since Xst is positive for both nucleons and A's, the total 5*°^ will also be positive 
for values of Tcut greater than 0.7 fm. That is, pion exchange between quarks produces a 
repulsive force within these baryons. 

To evaluate this contribution to the baryon energy, we must fix the value of gqqn-, as 



noted after Eq.([T]). To do this, we refer to results from the extended QDCSM [10(] and some 
general features relating the chiral coupling of pions to the axial current, /, to the coupling, 
(7, for the nominally equivalent (at low energies) pseudoscalar form (used here for its relative 
calculational simplicity). The usual relation at the nucleon level is 

QNN-k — In Nit l^-JOj 

where the value of fNNT,/{4:TT) = 0.075, as extracted from fits to nucleon-nucleon scattering 



14 



18| . The factor of 2Mi^ compensates for the fact that in the pseudoscalar coupling, 
ipj^ip, the upper and lower components of the Dirac wave function are coupled, with the 
lower being reduced, in momentum space, by a factor of order m^/ {E^ + Mn) ~ m^/ {2Mn)- 
The momentum scale is set by and the kinetic energy of the nucleon is small. 



In the extended QDCSM [10(], (which, as we noted previously, is a non-relativistic quark 
model of nucleon structure and nucleon-nucleon interactions extended from the QDCSM 
by the addition of pion exchange between quarks separated by a minimum distance), the 
relation between /nn-k and fqqT, is modified from the simple group-theoretic factor (5/3) by 
a small correction due to the quark delocalization within the nucleon, which accounts for 
the nucleon not being a point-like field lacking internal structure, viz., 

f 2 - A f 2 p-m2feV2 (Of.\ 
J qq-K — 2^-'^^^^ ■ ^ ' 

Here b ^ 0.6 fm describes the width of a single quark wave function in the nucleon and 
so provides for the experimentally determined root-mean-square size of the overall matter 
distribution in the nucleon. Since ^ 0.4, this correction is of order 10%, so that 

^ 0.025 . (37) 



Nonetheless, this value is not appropriate to our relativistic description here, as it has been 
extracted for non-relativistic wave functions. The effect is similar to the one evidenced in the 
right-hand-side of Eq. fl35|) . as our lower wave function components are suppressed relative to 
the upper components by a factor of order m^/(i?g + nig), the relevant ratio corresponding 
to what we saw above for the nucleon. Since the quarks here are highly relativistic, but with 
an energy of order one-third of the mass of the nucleon, the suppression is therefore not as 
great as for the case of the very non-relativistic nucleon. Since this suppression is already 
present in the scale of our lower wave function components, it follows that the appropriate 
value to use is 

g^jArr ^ 0.025 x fEl±^\" ^ o.lT , (38) 

where we have taken Eg = 361.6 MeV and rrig ~ 0, as in the LARQ. Note that this is 
approximately two orders-of-magnitude smaller than the value for nucleons, q^^^/Att, which 
suggests a more valid perturbation theory. 

The function /(rcut) reaches its minimum near r^nt = 1-05 fm, where its value is —0.492 
MeV. Multiplying together the various factors in Eq. ( !30l) . the maximum contribution for 
this exchange energy is 

Bl°l = 5.02 MeV, for the nucleon, 5*°*^ = 1.00 MeV, for the A. (39) 

This contribution to the self-energy produces a small change in mass, compared to rest 
masses of these baryons. It is comparable to electromagnetic corrections, which have been 
ignored at this level. 

This allows us to conclude that the fitted parameters of the LARQ model for the nucleon 
and A baryons are not significantly affected by inclusion of this pion-exchange contribu- 
tion between the quarks. The only remaining concerns are related to quark self-energy 
corrections, which we will deal with in following section. 

For larger values of r^ut, such as the 1.4 fm cutoff advocated in Ref. [l^, this contribution 
is even smaller. It is also interesting that, for a value of rcut ~ 0.7 fm (as suggested by the 
MIT bag model [19I), this contribution is exceptionally small, since /(rcut) crosses the x-axis 
near that value. 

However, none of this implies that the predictions of the LARQ model for nuclei, such 
as the deuteron, will be unchanged by this pion-exchange contribution. There, significant 
cancellations affect the net value of the binding energy. Hence it is still necessary to extend 



these calculations to the nuclear cases, and to determine the change in this pion-exchange 
interaction energy as the quark wave function changes in each particular nucleus. 



III. PION-LOOP CONTRIBUTION TO THE QUARK SELF-ENERGY 

In addition to the exchange of a pion between different quarks in the hadron, one must 
also consider, at the same level, the one-loop contribution to the quark's self-energy. (See 
Fig. [21) We first make an approximation that the major contribution to the loop integral 
for this self-energy comes when the intermediate quark state is given by the bound state 
used in the previous section. This approximation neglects the contributions of the excited 
bound-state quark wave functions. It is similar to one made by Inoue et al. jl^, but we 
improve on that work by keeping our calculation fully relativistic. 

We then compare the above bound-state approximation with two different approxima- 
tions, where we use a free massless quark propagator for the intermediate quark state. One 
version uses a sharp momentum cutoff, neglecting all pion momenta greater than kcut- The 
other uses a monopole form factor to cut off, more smoothly, pion momenta greater than its 
parameter, A. 

As we shall see, all three of these approximations give corrections on the order of a few 
MeV, i.e., are also small in comparison with the nucleon and A rest masses. 

Comparison of these three approximations weighs heavily on the question of whether or 
not this contribution can be calculated sufficiently reliably, so that the change that develops 
in a nucleus can be reasonably determined. 

A. Bound, Ground-state Approximation 

We begin by evaluating the first approximation to the loop correction. For a quark labeled 
"1" and a quark in the intermediate state labeled "2" , 




(40) 



where < r • r >= 3 and 



/ 



— m-jrf 



/""^(rcut) = E / d'n 



d^r2 {ipiil->ip2) 0(r - ^cut) (V'2^75V'i) • 



(41) 



Again, the B-function assures that r = |ri — > rcut . For the S'-wave quark wave functions 
of Eq. ([2]) we find 



iS2 



and 



(42) 



(43) 



Multiplying Eqs. (j42l) and (H3l) and using completeness of the spin function, J2s2 XS2XI2 — 1) 
we find 



51(^1^75^2) (^2^75^1) = 2 [^a(n)^b('^l) fi ■ f2 ^a('^2)^6('^2) 

+ |^a(ri)|V6(r2)r^ 



S2 



(44) 



where we have made use of the reality of the wave functions, the anticommutator {cxa, a/s} = 
6ai3, and the interchange symmetry/equivalence of ri and r2 to combine terms. 

With the pion Green's function factor, e~™'"''/r, the integral over the first term in Eq. 
( l44l) is two times /(rcut)? which we defined in the previous section and calculated to have 
the form given in Eqs. (JT^ and (l26l) . Thus, the self energy for a single quark from a pion 
loop consists of two terms. 



5££ = 6(^) [/(r,,0+/(rc 



(45) 



The new contribution to B^^, /(rcut), comes from the last term in Eq. (jSl). For a three- 
quark baryon, the total self-energy correction from pion loops will be three times B^^. 
We can evaluate / in a way very similar to that for the more complicated /. 



/(r, 



cut j 



-0(r - rcut) 



AttN k / x^dxi 0b(a;i) J(xi,a;cut) 
Jo 



(46) 



where 



J (xi , Xq 



d^X2 



-fix 



X 



-Q{X - Xcut) (I)1{X2) 



00 r c 

x^dx I dVtr — 



-fix 



X 



6(X - Xcut) 0a(|xi - X|) 



(47) 



Note that / is positive and there will be some cancellation in Eq. ( l45l) with /, which is 
negative for values of rcut of interest. How much cancellation depends on the relative sizes 
of / and /. 



Proceeding as in the previous section, with the Gaussian expansion of (pai^) defined in 
Eqs. (IMHD, we find, after some algebra, 

1 / IT V^"^ 

{e-^-erfc[titl(xcut)] - e+'^^^erfcftHlxeut)]} ■ (48) 

The quantities Cfc^fcj and ^^^^^(^cut) are the same as above, defined before Eq. (fT9!) and in 
Eq. fl25|) . Inserting J in Eq. fH6|) . we have finally 

/ \^^^ ~ 



fclfc2 

where 



^Ckik2 



IkM^cut) = Hxidxi <j)t\xi) {e-^"ierfc[4+i^(xcut)] - e+'''''erk[tlXM]} ■ (50) 

We have again evaluated this sum of one-dimensional integrals using Mathematica. It is 
a monotonically falling function ranging from 11.88 MeV at r^ut = 0.6 fm to 0.275 MeV at 
'"cut = 1-8 fm. For rcut = 1-0 fm, / = 6.17 MeV, quite a bit larger than / at this value of rcut 
(-0.487 MeV). Thus the sum of the integrals in the square bracket in Eq. (145|) is dominated 
by the / term. The solid curve in Fig. [3] shows the dependence of on r^ut- (The dashed 
curve shows the dominant contribution due to /(rcut)-) 

For Tcut = 1-4 fm that square bracket has a value of 1.436 MeV and the total value of the 
quark self-energy loops is B^°\i is 4.39 MeV. This is the contribution of all three quarks in a 
nucleon (or A) for this value of rcut- In comparison, for the nucleon, the exchange energy from 
pion exchange between the three quarks at the same value of rcut is, from Eq. fl30|) . 2.58 MeV. 
If this value were not much changed by the inclusion of higher energy intermediate states, 
it would be clear that the LARQ parameters are not significantly affected by the addition 
of "long-distance" pion exchanges. We next test this conclusion by using a free propagator 
approximation for the intermediate state. 



B. Free Propagator Approximation 



To compare with the above approximation we now evaluate, in four-momentum space, 
the pion-loop self-energy correction using a free-quark propagator in the intermediate state. 



This energy is, for a bound quark with four-momentum = {E, |p|), E = 0.362 GeV, 
an intermediate pion with four-momentum k^, and a massless intermediate free quark with 
four- momentum {p — k)^, 

d^p f d^k -J , l'^iP~k)u 



^i?.ee = ^9qq.j:j T^.J 7^.i^s{p)^l, 



X 



(27r)4y (27r)4""^^^^ [(^p ^ k) ^{p - k^ + te] 

{kf^k^^ — ml + ie) '^^'^^^■^^ 

= e/ MpHi-'i.i^siP) ■Mp.pn (51) 

The factor of 3 comes from the isospin factor, < r ■ r >. 
The loop integral Ji, defined in Eq. flSTl) is 

Mp,pn = / yfi u ,w" ^ • 1 (1.1., ^ 2 ^ • , = -^(P/^P") ' (52) 

J {27rp [{p - k) ^{p - k)i' + ie\ [k^k'' - + le) 
since p^, is the only available four- vector for defining the //-component. Its coefficient is 

^^""^ ' PuP" ' (27r)4 [p^p^^ - 2p^k^^ + k^k'' + le] {k^kf^ - ml + ie) ' ^ ' 

which will be evaluated below. 

With i^sip) = 4'l{p)lo and the factor of pu from J,^, the Dirac matrices in Eq. (ISTI) reduce 

to 

/ 

-Po cr • p 



7o75 YPu 75 



(54) 

ycr-p -po 

which is to be inserted between the spinors iplip) and ips{p)- We will see presently that this 
allows the spin sum over s to become trivial. 

So, there are two tasks: to determine the momentum-space wave functions ipsip) and to 
evaluate J{p^p^). 



1. The Momentum- Space Bound- State Wave Functions 

We need the four-momentum-space ips{p) corresponding to the bound-state Dirac wave 
function V'(x) given by Eqs. ([2]), (jl]), ([6]), and ([7]). To do this, we must add a time- dependence 
to 'V^(x), i.e., we need the four-space Fourier transform of ip{:x)e~'^^^ . 

The Fourier transform of the time dependence, e~*^*, is a bit tricky, as one needs to avoid 
the integration over the square of a 5-function, 5{po — E), in the integrals for -B|ee- However, 



normalizing instead in a box of time of size Tg, such that ue = ETe/{2tt) is an integer, we 
find the normahzed Fourier transform in time to be 



tPipOn) = {2t,TeY'^ 5n,nE . (55) 

Further, because of the discretization, integrals over dpo/{2TT) are replaced by sums involving 
the discrete eigenvalues pon = ^ttu/Te, 

r ^ V^nPo)/(po)^(po) - 7^ E i^\pon)f{27in/TE) ^{pon) = f{E) , (56) 

^-oo Zn 1e n=-oo 

The Fourier transform of the spatial wave function ^/'s(x) is more straightforward. Drop- 
ping the isospinor (which is already incorporated in the < r ■ r > factor in Eq. (15T|) above), 
we have from Eq. ([2]) 

V-(p)=/dVe-*-W = ivf2^V''( ^'^'l'^' I, (57) 



where g = |p|/k is dimensionless and 



k /ifc 

are also dimensionless real functions. 

With Eq. fl57|) . the wave-function dependence in Eq. fl5T]) becomes, using Eq. fl54|) . 

Y.^lipho75YPu7Mp) = -N' (^)' [po ^l{q) + 2|p| Mq)Mq)+po 4>l{q)] , (59) 

the spin-sum having become trivial, as advertised above. (The Pauli spin matrices again 
disappear from the problem since {<JiPi, CTjPj} = 25 ij PiPj = 2.) The three terms in Eq. (15^ 
are in decreasing order of importance. 



2. The Loop Integral 

The loop integral Jip^p^) is not the usual self-energy one-loop correction for a free particle 
for two reasons. First, the external quark is a bound-state particle, requiring an integration 
over its four-momentum, p^. The integration over dpo is trivial, using Eq. fl56|) . so the 



work to be done is in evaluating the triple integral over d\p\. (The four- momentum is 
time-like for |p| < E and space-hke otherwise.) Secondly, the integration over the pion 
four-momentum is "cut," corresponding to the above coordinate-space integrations being 
evaluated for r = |ri — > Tcut- Hence we want to restrict the integration over the pion 
three-momentum, |k|, to be from to kcut — ^c/rcut ~ 0.2 GeV/c. 

We first carry out the integral over dko (from — oo to -|-oo) using the Feynman trick, 

1 /"i 1 

-r^ dx 7 — — . (60) 

ab Jo [xa + {l-x)bf ^ ^ 

Letting a = (p — A;)^(p — k)'^ + ie, b — k^ki^ — rri^ + ie, and four- vector li, = ky — xp^,, we find 

J{p,Pn = I dx{l - X) ^ [y._i,2(^,,)p + ,, ■ (61) 
Here we have dropped terms linear in (they vanish) and have defined 

R''{p,x)^-x{l-x)p^p^ + {l-x)ml , (62) 

Note that R^, despite its appearance, is not positive definite; as a quadratic in x it can 
go negative for m^/x + < E'^. As we will see, this, among other things, leads to some 
imaginary parts in the evaluation of J{p^p^). 
Continuing, we have 

J{P,V^ ) = ldx{l-x) ^ y_ ^ (/g -S^ + ^ef ' ^^^^ 

where 

S\p,x,k) = \l\^ + R\p,x) . (64) 

Using Cauchy's integral theorem for integrating over the double pole at Iq — —S + ie, we 
find 

J(pX) -lf,dx{l- X) ^ [|,|2 + ^'(^^,)]3/. ■ (65) 

Unfortunately, the limits for the integration over d^l are angle dependent and complicated. 
Thus we go back to an integration over d^k. In the following we abbreviate the three-vector 
magnitudes |p| and |k| by p and k, respectively. 

J{E' - p2) ^ _1_ J^' dx (1 - x) K{p, x) , (66) 



where 



K{p, x) 

L{p,x,k) 
A{p, X, k) 
B{p, X, k) 



^cut 







hdk L{p, X, k) 



1 1 

-1 [A — B cos dy^i'^ 

2 I ™2 2 I r>2i 



2xpk > . 



(67) 

(68) 
(69) 
(70) 



Note that the functions A and B here are different from the Ak^^k2 and -Bfci,fc2 defined in Sec. 
lie More importantly, A{p^x, k) can be less than B{p,x, k) and even can become negative 
for certain small values of p and k when < a; < 1. Thus the integrals K{p, x) and L(p, x, k) 
can be complex in some regions of their variables. These imaginary parts arise in situations 
when the initial bound quark, sufficiently off-shell, can decay into an on-shell pion (with 
three-momentum k) and a free (massless) quark. Such decays will not contribute to the 
self-energy being calculated, since the quark is confined to the nucleon in reality. Therefore 
we discard any imaginary parts of K{p^ x) and L(p, x, k) below. 
Doing the integration over cos(^^), we obtain 



L{p,x,k) 



B 



1 



(71) 



VA^ VA + B_ 

where, as just discussed, it is understood that we only use the real part of L in the later 
integrations. In the limit that B ^ 0, L ^ 2/A^^'^, so there is no real singularity at /c = 0. 
The integral over dk defining K{p, x) can also be done analytically, yielding 



K{p,x) 
D{p,x) 

E{p,x) 
F±{p,x) 



D{p, x) + E{p, x) 



(kcut - pxy + R^ip, x) - J (fccut + pxY + R^ip, X 



2px 



F+{p,x)F_{p,x) 
R^ {p, x) 



2 



kcnt ±px + J {kcut ±px)^ + R^ip, x) 



(72) 
(73) 

(74) 
(75) 



In the limit that x — > 0, D{p, x) —2kcut/ y ^cut + ^'ii a-iid there is no singularity at x = 0. 
Also, noting that R^{p^l) = 0, one might be concerned about the logarithmic divergence 
at X = 1 of K{p,x). However, this singularity is moderated by the factor of (1 — x) in the 
integral over dx in Eq. (1661) . For the subsequent integration over p, the 1/p in D{p,x) will 
be cancelled by the p"^ in the d^p volume element. 



Again, we only take the real part of K{p, x) in the later integrations over x and p. For 
certain values of < x < 1 and (small) kcnt, the structure of K{p,x) can be quite complex. 
The argument of the logarithm in E{p,x) can pass through zero and even become complex 
in several places. Not only does the denominator R{p, x) become negative, but so can the 
arguments of the square roots in F±. An example for fccut = 0.100 GeV/c showing the real 
and imaginary parts of K is given in Fig. m 

It is also possible to do integrations over dx in Eq. flBBl) analytically, but the resulting 
expressions are extremely cumbersome and not very informative. Thus we have elected to 
do the remaining integration for J(£'^ ~ p"^) over dx numerically, along with the necessarily 
numerical integration over dp to obtain the value of Bf^^. 



3. Combining These Results 

From Eqs. (ISTI) . (156|) . (l59l) . (166|) . and (16711751) . we have, for a single quark, 

<e(fccut) = ^ I p'dp{E[<Pl{q)+<Pl{q)\+2pU<i)K<i)} 

X dx (1 - x)K{p,x) . (76) 

J 

As just stated, we have done the integrations over dp and dx numerically. We find that, 
as expected, there is a monotonic decrease in Bf^^{kcut) as the cutoff parameter fccut gets 
smaller (i.e., as rcut gets larger), as shown in Fig. O To make the comparison with Fig. [3] 
easier, we have plotted it as a function of Tic/ kcnt in fm. The fact that B becomes slightly 
negative about 1.4 fm is probably an artifact of our sharp cutoff approximation; see the next 
subsection for another version of the free propagator approximation in which the self-energy 
stays positive. 

If Tcut is identified with he/ kcnt (i-e., with proportionality factor 1.0), then the self-energy 
in this case is very small for values of rcut greater than 1.0 fm. The proportionality factor, 
however, need not be 1.0. 

For the nucleon (or the A), the result of Eq. fl76l) is to be multiplied by three. 



4- Alternative Loop Integral with a Monopole Form Factor Cutoff 



The loop integral J^, in Eq. fl52|) can also be "cut" by introducing a monopole cutoff 
factor —h?/ (A;^ — + it) in its integrand. This factor provides convergence to an otherwise 
logarithmcally divergent integral. It is normalized so its value is 1 when goes to 0. 
The momentum parameter A is (roughly) linearly related to the sharp cutoff /Ccut, and its 
proportionality factor will be discussed later. 

With the monopole factor in the loop integral above, JvipV)-, we rename it as 



d^k 



-A2 



(27r)4 [{p - fc) 2 + ^e ] (fc ^ _ ^2 + (^ 2 _ a2 + 



-A^ 



d^k 



A^-mlJ (27r)4 [(p - k) ^ + ie] 



1 



(A;^2 -ml + ie) {k^ - h? + ie) 



-A2 



A^ — rni 



■MPa,rnl)--MPu,^^) 



(77) 



Note that JuiP/^^^l) is just what we called Jvip^f) earlier. 
As before, with = rni and A^, 



J.(K^M2) = J(«^m2)p, . 



(78) 



The evaluation of the scalars J{p^,M^) proceeds much as before, but in fact the final 



integrations become simpler. With the same Feynman trick, 

dH 1 



(27r)4 [y^'-/^2(p^a;^M)]2 + ^e 



R'^{p,x,M) 



-x(l - x)p^ + (1 - x)M^ 



(79) 
(80) 



Integrating over the double pole in ko, we now find 



J{P',M' 



dH 



1 



"^""^^ (27r)3 [P + R^{p,x,M)]y^ 



(81) 



where / = Again, -)_ i?(p, x, M) can become negative for certain values of p and x. We 
continue ignoring any imaginary parts of Jlp^jM"^) when we calculate Ii^{p^). 

Doing the trivial angular integrations and integrating with respect to dl, we find the 
simple result that 



Up',^) 



A' 



IGtt^ yA^ — 
iPu 



K{p,A) . 



dx (1 — x) In 



- xp^ 

ml — xp^ _^ 



(82) 



where 



and we have defined 



A2 



rat 



dx (1 — x) In 



X 



X 



a{p,A) = A/p^' = A/{E'-p') , 
b{p,m^) = ml/p^^ = ml/{E^-p^). 



(83) 

(84) 
(85) 



[The function K{p,A) here is different from that defined by Eqs. fl66|) and fl72|) .] The 
logarithm in K has two singularities, when either A^—xp^ or rn^—xp^ equals zero. However, 
logarithmic singularities are integrable, so that is not a real problem. Assuming A > m^r, 
the argument of the logarithm in the integrand first goes negative when p"^ < — m^. As 
already discussed above, the resulting imaginary contributions to K{p, A) and hence luiPf^) 
are irrelevant for our self-energy calculation, and they will be discarded. 
Carrying out the integration over x in Eq. (l83l) . we find 

A^ r , . /a— 1\ / o \ . /a — 1 



Kip, A) 



A2 



a — 6 + In 



1 



2b] In 



- 2a] In 



(86) 



bearing in mind that a and b are functions of p. This result is valid even when the arguments 
of the logarithms go negative. For example. Fig. [H] shows the real and imaginary parts of 
K{p,A) as a function of p for a value of A = 0.350 GeV/c. Note the smoother behavior 
of the real part as p goes through the singularity at p = J E'^ — m^, where /^.'s value is 



ln(A/m^). For smaller values of A the real part is more sharply peaked at E"^ — and 
the width of the positive imaginary "spike" between the two singularities is narrower. 
From Eqs. (|HTi) . fISB]) . (13^ . and (15^ . we have, for a single quark. 



Bll{A) 



3 9qq-K 



p'dp{E [^l{q) + 4>l{q)]+2pMq)Mq)}K{p,A) . (87) 



2 Air JO 

We have carried out the integration over dp numerically. There is, in this approximation, 
again a monotonically decreasing dependence on the inverse of the cutoff parameter A, as 
shown in Fig. [71 (Here we have chosen x-axis variable to be he/ A, in units of fm, to compare 
better with the ground state approximation shown in Fig. [3l) 

The comparison of Fig. [7] with Fig. [5] was at first surprising to us, but the results for the 
value of are comparable if A ^ 2.5/i;cut, as shown in Fig. [HI 



C. Comment on renormalization 

The reader familiar with field theory may be concerned that we have not carried out mass 
and wave function renormalization subtractions of the quantity calculated above. We have 
not done so because the renormalized field theoretic degrees of freedom in the LARQ are 
massless, current quarks. To the extent that pion self energy corrections must be applied, 
we view these as implementing the construction of constituent quarks, which have revised 
properties, as remarked in the Introduction. 

Nevertheless, the fundamental approach of the LARQ was to carry out calculations with 
current quarks and with minimal modeling of QCD. From that point of view, the pion loop 
self-energy correction should be subtracted and the quark wavefunction correspondingly 
renormalized to recover the current quarks of the model. If this is carried out then, since 
the calculations above produce relatively small values, this makes clear that the changes after 
renormalization must be very small. We may then further infer that the changes induced 
by the nuclear distortions of the nucleon quark wavefunctions produce differences between 
a nucleus and A nucleons that are even smaller. 

It thus appears to be reasonable to conclude rather strongly that, as for the pion exchanges 
between quarks analyzed in SecUll these effects on the parameters of the LARQ due to 
including pion-quark coupling are negligible. The only question remaining for the LARQ 
is how much the pion interactions between quarks from different nucleons affects the total 
binding energy of the nucleus. We will address this question for ^He and ^He in a later 
publication. 

IV. RELATION TO OTHER CALCULATIONS AND CONCLUSIONS 

The calculations in the previous two sections correspond to calculation of the two graphs 
in Fig. m The exact result for Fig. [2)d requires a summation over all intermediate excited 
states of the quark within the confining potential or, equivalently, use of the in-medium 
quark propagator rather than the free propagator. However, as the two results of using only 
the lowest bound state and the free propagator are both small, we conclude that we have 
a reasonable estimate of the long-distance pionic contributions to the mass of the nucleon 
and of the Delta within the LARQ quark model of these states. 



Estimates of the pionic (pion cloud) contributions to these states have also been made 
using chiral pertubation theory (xPT). These calculations are restricted (so far) to the 
contributions shown in Fig. [2] in a form which ignores the substructure and uses the low 
energy constants (LEG) that are extracted from data using xPT. It should be noted that 
the LARQ calculation also implicitly includes all other intermediate states in the pion loop 
calculations of xPT - all A^*, A*, etc. states. From xPT, one expects their contributions 
to be small and decreasing with increasing mass. As such, it is interesting to compare xPT 
results with those of the LARQ. 

Unfortunately, there does not seem to be much agreement among the various xPT cal- 
culations as to the result. For example, Ref. ^] obtains a large correction 200 MeV) 
for the total nucleon self-energy, using only the ground state baryon as the intermediate 
state, (where the substructure is used to compute an effective vertex function for the nNN 
coupling rather than working with the substructure directly as done here). Other papers 21 1 
include A and A^* intermediate state resonances with ambiguous results. 

In conclusion, we find that the contribution of this additional interaction to the mass 
of the nucleon (and the A-baryon) is small. The scale of effects that we have found is 
comparable to the expected scale for electromagnetic corrections, which were not included 
in the LARQ and so contribute to its overall scale uncertainty. We conclude from this that 
the LARQ is stable under this change, in the sense that significant changes to the model 
parameters are not required to maintain a match to data on state masses. We also infer 
that inclusion of these pion exchange contributions in calculations of nuclear states will not 
significantly affect the quark structure of the nuclei considered, although it could improve 
the agreement with actual binding energies. We will report results of those calculations 
elsewhere. 

This work has been supported by the Department of Energy under contract W-7405- 
ENG-36. We acknowledge useful conversations and correspondence with J. L. Friar, J. N. 
Ginocchio, A. K. Kerman, M. M. Nieto, Jialun Ping, K. E. Schmidt, and A. Steiner. 



V. APPENDIX: GROUP THEORETIC CALCULATION OF THE SPIN-ISOPIN 
SUM 



An elegant (and relatively quick) calculation of the Xst factor uses some group theory 
techniques. First, rewrite Xst in terms of the spin and isospin angular momentum operators, 

Xst = 8 ^ < 5, M5I s^*) ■ s^^'^ \S, Ms >< T, Mt\ t^'^ ■ t'-^^ |T, Mt > . (88) 

The reason for doing this is because the three quarks are in a symmetric representation of the 
group SU(4), for which the generators are Sa, ta, and Sat/3 in the fundamental representation. 
For our case of three symmetric quarks, these generators become 

33 3 



i=l 1=1 i=l 



S^ = Y: , = E t'a^ , and G^p = ^ s^kf . (89) 



The sum in Xst can be evaluated from the quadratic Casimir invariant €2^^ for the 
symmetric SU(4) representation for and A, whose Young tableau has "overhang numbers" 
(A,jU, z/) = (3,0,0). For a general SU(4) representation 



22 



C2^^ = ^[3A(A + 4)+4/i(^ + 4) + 3z/(z/ + 4)+4^(A + z/) + 2Az/] (90) 
= 63/4, for (A,/x,z/) = (3,0,0) . (91) 

The C^') can also be expressed in terms of the total spin and isospin quantum number and 
the scalar formed from the Gap as [23,] 

= S{S + 1) + T(T + 1) + 4 < GapGap > , (92) 

Breaking up product in the last term into terms with i = j and those with i ^ j, 

G ^G ^ - S^Sf^^Sf^^ + V 

^aP^aP — ^S^lp S^lp + 2^ Ip 
i=l i¥'j=i 

= f](s«.s«)(t»-t») + lx5T. (93) 

i=i ^ 

The sum on i in the first term is easy. 



Thus 



X 



27 



ST 



< cj'^ > -s{s + 1) - r(T + 1) 



18-2[5(5 + l)+T(T + l)] 



using the value of 63/4 for < 6*2^'' > given in Eq. (iQli) . The spin-isospin sums for and A 
are therefore 

Xl/2, 1/2 = 15 , X3/2, 3/2 = 3, (96) 

in agreement with the values given in Sec. IIIDi 
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FIG. 1: Variation of integral for exchange between different quarks, /, as a function of rcut- The 
X- and y-axes have units of fm and MeV, respectively. 
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FIG. 2: Both the self-energy pion-loop, and pion exchange between separate quarks, appear as 
self-energy corrections to the nucleon or the A when the internal substructure of the baryons is 
ignored, as in Chiral perturbation Theory. 




FIG. 3: The solid curve shows the pion-loop self energy correction for a single quark as function 
of Tcut- The dashed curve shows the dominant contribution coming from I{rcut), the difference 
coming from the smaller /(rcut) displayed in Fig. [TJ 
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FIG. 4: K{p,x) for x = 0.3 and A^cut = 0.100 MeV/c as a function oi p (in MeV/c). The solid 
curve displays the real part of K, the dashed curve is for the imaginary part. The "spike" in the 
imaginary part near p = 0.26 MeV/c is actually a narrow flat-topped pedestal with value +7r. 




FIG. 5: The free propagator approximation for the self-energy of a single quark with a pion-loop 
correction with sharp cutoff, fccut- ^ft.ce(^cut) is shown as a function of hc/kcut for comparison with 
Fig. El 




FIG. 6: The real (solid curve) and imaginary(dashed curve) parts oi K{p, A) for a value of A = 0.350 
GeV/c, close to, but less than, the bound quark energy, E = 0.362 GeV. 




0.6 0.8 1 1.2 1.4 1.6 1.8 

nc/A (fm) 



FIG. 7: The self-energy contribution for a single quark, i?|,gg(A), for a pion-quark loop correction 
using a monopole cutoff parameter A. It is shown as a function of Tic/ K in fm, for comparison with 
both Fig. [3 and Fig. [SI 
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FIG. 8: The soUd curve shows i?|,gg(A), as in Fig. [71 while the dotted curve is B^^^ ( ^cut ) ? in 
Fig. El The dashed curve indicates how the agreement between these two curves for A reduced by 
a factor of 2.5. 




